Distribution of proper delay times in quantum chaotic scattering: 
A crossover from ideal to weak coupling 



O 

o 

(N 

< 
(N 



Hans-Jiirgen Sommers 1 , Dmitry V. Savin 1,2 , and Valentin V. Sokolov 2 
1 Fachbereich Physik, Universitat-GH Essen, 45117 Essen, Germany 
2 Budker Institute of Nuclear Physics, 630090 Novosibirsk, Russia 
(Received 8 May 2001; published 8 August 2001: Phys. Rev. Lett. 87, 094101 (2001)) 

The probability distribution of the proper delay times during scattering on a chaotic system is 
derived in the framework of the random matrix approach and the supersymmetry method. The 
result obtained is valid for an arbitrary number of scattering channels as well as arbitrary coupling 
to the energy continuum. The case of statistically equivalent channels is studied in detail. In 
particular, the semiclassical limit of infinite number of weak channels is paid appreciable attention. 
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The temporal aspect of collision and transport phe- 
nomena has, starting with the pioneering works re- 
peatedly attracted attention. Analysis of the duration 
of a process gives an interesting and important informa- 
tion complementary to that delivered by the energy rep- 
resentation. The temporal approach becomes especially 
elucidating in the cases of chaotic resonance scattering 
and transport in disordered media when many compli- 
cated long-lived intermediate states are involved [|-@; 
see for recent reviews. 

The delay of an almost monochromatic wave packet 
during an M-channel process is conventionally described 
by the MxM Wigner-Smith matrix Q^-ihS^dS/dE, 
with S being the scattering matrix at the collision energy 
E. In the resonance scattering, the matrix element Q cc i 
describes the overlap of the internal parts of the scatter- 
ing wave functions in the incident channels c and c' [ p"8| . 
In particular, the diagonal element Q cc , which is inter- 
preted as the mean time delay in the channel c, gives the 
norm of the internal part b c (E) of the scattering wave 
function excited through the channel c (the normaliza- 
tion is fixed as usual by setting the flow in any given en- 
trance channel to one). This directly relates the Wigner- 
Smith matrix to the effective non-Hermitian Hamiltonian 
7i=H—{i/2)VV^ of the unstable intermediate system as 
follows (henceforth Ti=\) fl8fl: 



Q(E) = 
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V = b\E)b(E) 



The Hermitian part H determines the appropriate basis 
for the internal motion whilst the amplitudes V£ describe 
the coupling between N interior and M channel states. 
We adopt below the standard statistical approach |p| 
based on the random-matrix theory (RMT) to simulate 
the complicated intrinsic motion; see 20 21J. 



Recently, an appreciable success has been achieved by 
the authors of Refs. |l2|. By quite a distinguished way of 
reasoning, they managed to calculate in the framework 
of RMT the joint probability distribution of the eigen- 
values q c of the Wigner-Smith matrix, which are called 



the proper delay times. Using then the method of or- 
thogonal polynomials they calculated also the density of 
proper delay times. These findings received later experi- 
mentally relevant applications to quantum dots j2i| and 
optics of random media . However, all the results are 
restricted to the case of perfect coupling to the contin- 
uum (all transmission coefficients are equal to 1). 

In this Letter we show that representation (|l|) allows 
us to extend at least part of the results mentioned above 
to the case of arbitrary coupling. Namely, we calculate 
the probability distribution of the proper delay times q c 
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where (. . .) denotes the statistical average, and 

G(z) = M- 1 (tr {z - Q)- 1 ) (3) 

is the trace of the Green's function of the Wigner-Smith 
matrix. The positive definite form ([!]) of Q provides the 
analyticity of G(z) throughout the complex z-plane save 
the discontinuity line along the positive part of the real 
axis q=Ke z>0. Function (^) is normalized to unity. 

The function G(z) can be readily obtained from the 
generating function (z±—z±u>/2) 



(1) Z(z,w) = 
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as G(z) = M~ 1 dZ(z, u>)/duj\ u —o. We have exploited here 
the explicit form (|l|) to cast Z(z, ui) in terms of determi- 
nants of Nx N matrices A± (rather than those of MxM 
ones z± — Q). Now, we use the identity 
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and double the dimension, introducing Pauli matrices o~i 
in the doubling space. As a result, the determinants get 
the following simple form: 



1 



Det(At) - Det(iw t - ^-(l + oi) - i(f7-ff) 03), (6) 

which allows us to represent the generating function (Q) 
in form of a Gaussian integral over an auxiliary super- 
vector and therefore to e mploy hereafter the standard 
supersymmetry technique [[l9| , |24| . As a result, we come 
to the following (exact at N — > 00) expression for the 
generating function 

Z( z , ,)= -± / [da] cxp [— str(*7(A - iA x )a) 

M 

x JJ sdet[l + j c A(iE/2 + Tiv(E)a)]- 1/2 (7) 



as the integral over a noncompact saddle-point manifold 
<T. The 8x8 supermatrices A, Ai appearing above are 
the supermatrix analog of the Pauli matrices 03, u\, and 
L f =diag(zZ 1 l2, z+ 1 l2, -zZ 1 l2, 2;+ 1 l2)- Definitions of the 
superalgebra as well as the explicit parameterization of a 
and the invariant measure can be foun d in |0|§g . The 
average density of states v{E)=^ x ^J\— (E/2) 2 deter- 
mines the mean level spacing A-(iyN)^ 1 of the closed 
system. Phenomenological constants 7 C >0 originate from 
the coupling amplitudes: (V^V) cc i—2^ c 8 cc i\ they enter 
final expressions only by means of the transmission co- 
efficients T c =l- |5 cc | 2 =2[1 + ( 7c + 7 c - 1 )/2t W (£;)]- 1 < 1 
|^9| . The latter are the conventional characteristics of 
coupling of an unstable system to the energy continuum. 

Suspending all technical details of quite standard cal- 
culations [for the sake of simplicity, we restricted them 
to the unitary symmetry class which corresponds to sys- 
tems with broken time-reversal symmetry (TRS)] to a 
more extended publication, we proceed with the follow- 
ing result for the Green function: 
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Henceforth, we scale the original variable z—C^tn [and 
correspondingly the Green's function (^)] in the natural 
units of the Heisenberg time tH=2irfr/ A. The constants 
g c = 2/T c — 1 are related to the transmission coefficients 
T c . At last, we denote &(Ai)=e Al / c I (C~ V A i = 1) and 
/(A 2 )=e- A2 /«J (C- i yT^i), with l (x) (Mx)) being 
the modified (usual) Bcsscl function of zero'th order. As 
it stands, the Green's function (^|) is an analytic function 
of the complex variable £ in the half-plane Re£<0. One 
can deform the original Ai-integration to that along the 
imaginary axis [1,+zoo) [or [1,-ioo)] for Im(C _1 )>0 [or 
lm(£ _1 )<0], so that G(C) is analytically continued to the 
whole complex £-plane with a cut along positive Re 
Applying (Q), we get finally the following expression for 
the probability distribution of the proper delay times 
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for t>0, and VM(t) = for negative times. Here 
F C B = e->»' t Ut- l y/gp^) [] . (10) 



9a~ 9c 



F? = U + dXe-V* Mt-'Vl^) II (-9« + A )- ( n ) 

The function (O) is, in fact, a polynomial in 1/t. The 
formulae (P)-(pT[) are valid for an arbitrary number of 
scattering channels M and arbitrary constants g c > 1. In 
the situation when all chaotic states of the target sys- 
tem are alike, the channels get statistically equivalent. 
Setting all g c =g, we arrive in such a case at 
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with 



Bm = ±(-^) M e-3/n (t- l ^/g^i), (13) 
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It is worth noting that, due to similar analytical structure 
of the functions Ff and Ff , distribution (||) decays as 
£— M-2 m limit >oo. This is in agreement with 
the known universal law t~W2-a gg^-g^j valid for all 
Dyson's symmetry classes (3=1,2,4 ((3=2 in the present 
case of broken TRS). 

In the single-channel case, M=l, Eqs. (||) and (JlJ) 
confirm the result Vi(t)=t~ 1 dBa/dt found earlier in pfjj . 
When the constant g 2 /T 3> 1 this distribution shows 
a narrow maximum at the small time t ~ (2g) _1 <C 1 
because of the cooperative influence of the tails of many 
remote resonances. The far asymptotics 'Pi(i)cxt~ 3 , t^> 
g, is due to fluctuations of the widths of many narrow 
resonances HUTU- At last, in the parametrically wide 
domain (2g)^ 1 <^.t<^g the distribution decreases as t~ 3 / 2 
M , Q| . Just this domain saturates in the case g ^> 1 
the well known [[| sum rule (t) ={trQ) jtn=\, which is, 
actually, satisfied regardless of the coupling strength to 
the continuum [see Eq. (|l^) below]. 

At the same time, the average time spent by a spa- 
tially small wave packet in the interaction region (colli- 
sion time) is equal to (r) = (Q)/T = t H (t)/T §||f|[l7|. 
This relation suggests splitting the transmission coeffi- 
cient off the delay time of an almost monochromatic 
wave: q=Tr. The new variable r is then the time the pro- 
jectile is delayed inside the interaction region after having 



2 



penetrated there with the probability T . In agreement 
with the above consideration, the most probable dura- 
tion r of the stay inside is given by the Heisenberg time 
t H . 

In the case of an arbitrary number of channels M an 
additional scaling t s =Mt provides a similar condition 



(ts 



POO 

qM/t H =Tt)= / dt s t s V s (t s ) = 1 . 
Jo 



(15) 



Here V s (t s )= M _1, Pftf (i=t s /M) and the mean transmis- 
sion coefficient T=M~ 1 ^2 c T c has been set off, q=Tr. 
The scaled variable t s measures the time r in units of 
a typical life time. Indeed, the characteristic width of a 
state which decays through M open channels is estimated 
asr=£T c A/27r 0. 




t s -Mt 

FIG. 1. The distribution Vs(t 3 ) of the scaled proper delay 
times at M=2 and T=0.25 for the case of equivalent, ST=0, 
and strongly nonequivalent channels, <5T=0.49, (dashed and 
solid lines, respectively). 

Figure 1 demonstrates, using the example of the 
two-channel distribution V%{^)^ the dependence on the 
strengths of individual channels. To fetch out the in- 
fluence of channel nonequi valence, we keep the openness 
characterized by Ti+T2=2T fixed and change the differ- 
ence \T\— T2\=5T . Two distinct time scales with similar 
probabilities are well seen when ST is close enough to 
its maximal value min(2T,2— 2T). Fast processes appear 
when one of the transmission coefficients becomes small. 

We consider now in detail the case of large number of 
statistically equivalent channels, M^>1. The regime of 
isolated resonances, T^;A, corresponds to the condition 
g^> M>1. In this limit the widths of resonances become 
statistically independent and expression ( [l2|) reduces to 
the probability distribution of the partial delay times in- 
vestigated in much detail in (l(J , see also Jll],|l6| . 

Here we focus our attention rather on the scmiclassical 
limit of strongly overlapping resonances excited through 
and decaying into a number M of channels, which is 
scaled with N, the strength constant g being kept arbi- 
trary. Similar to Refs. ^5j,|), an additional saddle-point 
approximation can be used in this case to perform the 



integration in (Q). Under the physically justified con- 
dition H M/iV<l (though both M, iV— >oo) one finds 
that P(t„) = WmM^ooVsits) = {irtl)~ l 1xQ.K{t a ), where 
the function K(t s ) = t 2 (G(t s )—tJ 1 ) satisfies a cubic equa- 
tion of the form 



K(K 2 + 2gK+ 1) -t s (K 2 - 1) = 0. 



(16) 



Being derived for the case of the unitary ensemble this re- 
sult remains actually valid for all tree symmetry classes. 
The calculation just described sets a new intrinsic time 
scale of the problem j2j| : the empty gap r g between the 
real axis in the complex energy plane and the cloud of 
the resonances in its lower half plane. In the limit consid- 
ered T g coincides exactly with the well-known Weisskopf 
width T w =MTA/2ir. 

The case of ideal coupling, g=l, is especially sim- 
ple since the cubic equation (|l6| ) readily reduces to a 
quadratic one which immediately leads to the result 



P id (t s ) = (2irt 2 s )-\ (t*-t s ){t 



(17) 



with ii^=3±v / 8, obtained earlier in jl2| from the Laguerre 
ensemble. For an arbitrary value of g the searched dis- 
tribution is expressed as 
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(18) 



in terms of two polynomials: D=— 3r/' i + 2grf + llr) 2 /3 — 
4.9?7 + .9 2 + 3~ 3 and r=T] 3 ~2rj + g, where -q = (2g — t s )/3. 
The density does not vanish in the domain where the 
discriminant D > 0. It is obvious that the forth order 
polynomial D always has two real roots r\- < and ?7+ > 
and D is positive only inside the region r\- < rj < 77+. 
Therefore, similar to the case of perfect coupling, the 
time t s is always restricted to a finite domain t- <t s <t+. 
Near the edges t±=2g ±3\rizf\ the discriminant D is small 
and expression (|l8|), as in (|l7|), simplifies to 



P{t B ) « {■Kt 2 s )- 1 \r\- i l*^D/Z<x y/\t^t^ 



(19) 



A more detailed investigation is possible in the case of 
week individual channels, g^>l- All the roots of D in this 
case can be found as expansion due to the smallness of 
1/g. Near the lower edge t-~l/8g Eq. (pi) gives 



P(ts) 
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(20) 



In this case the term |r| dominates in Eq. fllSj ) and there- 
fore Eq. ( p0| ) is valid in the parametrically large do- 
main I/85 < t s 2g— |(5/2) 1 / 3 where the influence of 
other (three) roots remains negligibly weak ^(|. When 
t s I/85 this distribution displays the universal behav- 



ior t 



-3/2 



Such a law is a typical and most robust feature 



of the time-delay distributions in a weakly open chaotic 
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system and may be expected from quite a general argu- 
mentation |lO|JTl| ]. Within a relatively narrow domain 
-K5/2) 1 / 3 <i s -2g< 3(3/2) 1 / 3 the distribution devel- 
ops a local minimum and maximum, which are due to a 
pair of the complex roots i ± «2g-|(.g/2) 1 / 3 (l±iV3). At 
last, very close to the upper edge t + ~2g + 3(g/2) 1 / 3 the 
distribution sharply changes to the form 



(21) 



V3irt 2 s 

provided by Eq. (|l9|). Figure 2 illustrates these results 
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FIG. 2. The distribution P(t a ) in the limit M-^oo. The 
exact ( |l^ ) and approximate expressions @ (dotted) and (§]]) 
(dashed line, see inset) are shown at p=10. An intermediate 



t. 



-3/2 



behavior is clearly seen. 



With the help of the distribution found the Green's 
function can be represented in the spectral form G(z) = 
L + dt s P(t s )/(z ~ t s ). This function is analytical both 
at the origin and infinity and therefore can be expanded 
near these points in power series. Such an expansion 
yields the general formulae for the negative and positive 
moments of the distribution (n>0) 



(n+l)\_ ~1 d n G(z) 



dz n 



1 g Gfr 1 ) 



C=o 



All the moments are finite. The derivatives can directly 
be calculated by making use of the cubic equation (|l6|) . 

In conclusion, we have derived the distribution of 
proper delay times in chaotic scattering at arbitrary num- 
ber of scattering channels and arbitrary coupling to con- 
tinua. The appearance of distinct time scales is traced 
to the properties of the delay time distribution with in- 
creasing degree of channel nonequi valence. The physi- 
cally interesting case of many equivalent channels weakly 
coupled to continua is studied in much detail. 
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